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The branching ratio of the electromagnetic rare decays 77 — > %®yy and tj' — > (ft , 77)77 are ana- 
lysed in terms of scalar and vector meson exchange contributions using the frameworks of the 
Linear Sigma Model and Vector Meson Dominance, respectively. The measured 77 — > n°yy pro- 
cess serves as a test of our approach while the non yet measured 77' — > (ft , 77)77 reactions are 
predicted for the first time. Our prediction for the 77 — > n°yy decay agrees with recent experimen- 
tal reported values, thus supporting the validity of our framework. Therefore, our predictions for 
the 77' — » 7T°yy and 77' — > 77 77 decays should be taken as a first indication of the possible values of 
the associated branching ratios. We hope these predictions to be interesting and useful for experi- 
ments such as KLOE-2, Crystal Ball, WASA, and BES-III where these processes are expected to 
be measured in the next future. 
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1. Introduction 

The electromagnetic rare decays rj ' — > n°yy and 77 ' — » rj yy are calculated for the first time 
with the aim of completing existing calculations on the related r\ — > n°yy process. The first two, 
which have not yet been observed, will be potentially measured at several experiments such as 
KLOE-2, Crystal Ball, WAS A, and BES-III [Q, @|, while the latter, for which a measurement of its 
branching ratio and invariant mass spectrum already exists will be certainly determined with 
higher precision. 

From the experimental point of view, the situation at present is the following. The branch- 
ing ratio (BR) of rj ->• n°yy has been measured by GAMS-2000 [§, BR = (7.2 ± 1.4) x 1(T 4 , 
and CrystalBall@AGS in 2005 [g], BR = (3.5 ±0.7 ±0.6) x 10~ 4 , and 2008 [g], BR = (2.21 ± 
0.24 ±0.47) x 10~ 4 , the latter also including an invariant-mass spectrum for the two photons. The 
PDG 2010 fit is BR = (2.7 ±0.5) x 10~ 4 [§]. More recently, preliminary results from Crystal- 
Ball@MAMI |0|, BR = (2.25 ± 0.46 ± 0.17) x 10~ 4 , and KLOE |§ g|, BR = (0.84 ± 0.27 ± 
0.14) x 10~ 4 , have been reported as well. For the tj' — > n yy decay, only an upper bound exists, 
BR < 8 x 10~ 4 at 90% CL, obtained by the GAMS-2000 experiment 0] 25 years ago. Finally, for 
rj' — > rj yy there is no experimental evidence so far. On the theory side, the T7 — > n°yy process has 



been studied in many different frameworks, Chiral Perturbation Theory (ChPT) Ql 1Q, constituent 



quark model [12], three-flavor Nambu-Jona-Lasinio model [|||], a chiral unitary approach [|l4|, |l5| ], 
most of them in combination with the Vector Meson Dominance (VMD) prediction. On the con- 
trary, there are no theoretical analyses neither for t}' — > 7T°yy nor for rj' — > rjyy. 

It is the purpose here to give an estimate of the branching ratio of these three processes. Since 
we are more interested in an estimate rather than a detailed calculation, we will include in our 
analysis only the two main contributions, that is, the exchange of an intermediate vector meson 
through the decay chain P° — > Vy — > P°yy plus the chiral loops. Later, the chiral-loop prediction 
will be substituted by a Linear Sigma Model (LaM) calculation where the effects of scalar meson 
resonances are taken into account explicitly. As a check of our approach, we first calculate these 
two contributions for the case of T7 — > 7i°yy. Then, for the first time, we perform the same analysis 
for t]' — > 7r°yy and r\' — > r\yy. 



2. Chiral-loop prediction 

In ChPT, the tree-level contributions at &{p 2 ) and &{p A ) vanish because the pseudoscalar 
mesons involved are neutral. The first non- vanishing contribution to 77 — 5- 71° yy comes at 0(p 4 ), 
either from loops involving kaons, largely suppressed due to the kaon masses, or from pion loops, 
again suppressed since they violate G parity and are thus proportional to m u — mj. Numerically, it 
is seen to be three times smaller JTT|]. To simplify, we neglect the second contribution and work in 
the isospin limit. The first sizable contribution comes at \ but the coefficients involved are 
not well determined and one must resort to phenomenological models to fix them. In this sense, 
for instance, VMD has been used to determine these coefficients by expanding the vector meson 
propagators and retain the lowest term. This leads to values for the T7 — > 71° yy decay rate two 
times smaller than the "all order" estimate keeping the full vector meson propagator JTT|] . For the 
same process, the contributions of the scalar ao(980) and tensor «2(1320) resonances to the ^(p 6 ) 
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chiral coefficients were calculated in the same manner but no "all order" estimate was given in any 
case. In addition, contrary to the VMD contribution where the coupling constant appears squared, 
the signs of the «o and ci2 contributions are not unambiguously fixed. On general grounds, one 
would expect the 02 effects to be smaller than the ao ones due to the heavier mass involved in the 
propagators. For this reason, we will consider only the scalar meson contributions to the processes 
under analysis and provide an "all order" estimate of these scalar effects based on a calculation 
performed in the LaM model. In this way, we will be able, first, to fix the sign ambiguity and, 
second, to test the relevance of including the full scalar meson propagators, in a given model, 
instead of integrating them out. However, for the sake of completeness, we start considering the 
dominant chiral-loop contribution, that is, the contributions containing two vertices of the lowest 
order Lagrangian and a charged pion or kaon loop. The <^(p 8 ) loop corrections from diagrams with 
two anomalous vertices are seen to be very small JTT| ] and thus not considered here. The explicit 
contributions of intermediate vector and scalar mesons are postponed to the next sections. 

We start discussing the 17 — > 7T°yy case. As stated before, the contribution from kaon loops is 
dominant and the pion loops vanish in the isospin limit. The amplitude is written as 



2a 1 



71 m 



'k- 



(2.1) 



where {a} = (si -£2){q\ ^2) — (£1 ^2) (£2 'li), £1,2 and are the polarization and four-momentum 
vectors of the final photons, sk = s/m 2 K+ , s = {q\+ q^) 2 = 1q\ ■ q2 is the invariant mass of the two 
photons, L(s) is the loop integral defined as 
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and ^z+jr-^jfl is the four-pseudoscalar amplitude 
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with (pp the 77-77' mixing angle in the quark-flavour basis, resulting from the loop computation 
(not to confuse it with the four-pseudoscalar scattering amplitude calculated in ChPT at lowest 



order). It is important to notice that in the seminal work of Ref. [11] this chiral-loop prediction 



was computed taking into account the rj & contribution alone and the mixing angle was fixed to 
dp = (pp — arctan\/2 = arcsin(— 1/3) ~ —19.5°. Now, the T70 contribution is also considered (in 
the large-A^ c limit where the pseudoscalar singlet is the ninth pseudo-Goldstone boson) and the 
dependence on the mixing angle is made explicit. 

For the r\' — > 7T°yy case, the associated amplitude is that of Eq. ( |2.1[ ) but replacing — > m^, 
(cos (pp + y/2sin(pp) — > (sin(pp — a/2 cos (pp) and (cos (pp — sin<pp/\/2) — > (sin (pp + cos (pp/ \/2) in 
Eq. @. 
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Finally, for the v\' — > r\yy case, two amplitudes contribute, one through a loop of charged 
kaons, as in the former two cases, and the other through a loop of charged pions, which in this case 
is not suppressed by G-parity. Again, the corresponding amplitudes are that of Eq. (^~l|), replacing 
sk — > s K and m K + — > m K + for the pion loop, with 



K+K-^r\r]' 



1 

"4/2 



+-(2mj c — m 2 c ) ( 2sin2<pp 




r- _ sin2<pp 
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and 



2^sin2<p P 



(2.4) 



(2.5) 



The latter amplitude coincides with that of r\' — > 77 7T+7T when computed in the large-rV c ChPT at 
lowest order [16]. Needless to say, the former amplitudes for 77' — ^ n° yy and 77 ' — >• 17 77 constitute 



the first chiral-loop predictions of these two processes. 



3. VMD prediction 



Next to the chiral-loop amplitudes, there are also "all order" estimates of the corresponding 
exchange of intermediate vector bosons which are calculated in the framework of VMD. The full 
VMD amplitude was seen to produce the dominant contribution to 77 — > 71° 77 [11], and the same 
happens, as we see below, for r\' — > n°yy and r\' — >■ 77 yy. Now, we review the calculation for the 
77 — > 71° yy case, with some improvements with respect to Ref. [jTiQ, and then calculate for the first 
time the full VMD amplitudes of t]' 
as 



71° 77 and r\' — > 77 yy. For 77 — > n Q yy, the amplitude is written 
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where t,u= {P — q 2 ^Y 



■ mt 



-2P-q 2A ,{b} = (e l -q 2 )(e 2 -P)(P-q l ) + (e 2 -q l )(e l -P)(P-q 2 )-( y e v 



e 2 ){P-qx){P-q 2 ) - (£1 ■P){e 2 -P){q l -q 2 ) mxdD v (t) 



mt 



gators for V = ft), <f) . For the p we use instead an energy-dependent T p (t) = T p x [(? — 4m|) / (m p 



imyTy are the vector meson propa- 

Yp ^ a™- 2 ^ it™- 2 

4m 2 )] 3 / 2 x 6(t — 4m 2 ). For r\' — > 7C°77and 77' — > r\yy, the related amplitudes are Eq. ( |3.1| ) with the 
replacements gv^ygvifiy ~* gv-q'ygvTpy and 8vr)ygvrf>y -> gv-q'ygv-qy, respectively, and m 2 -> m 2 ,. 
When the OZI-rule is applied, that is ft) = (uu + dd)/^/2 and = sj, the corresponding couplings 
are 

0, 



gprj'ygpTiPy = gwq'ygmjfiy : 
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where G = 3g 2 / (4it 2 f n ) and g is the vector-pseudoscalar-pseudoscalar coupling constant of VMD 
which can be fixed from various p and ft) decay data. 
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chiral loops 


LaM 


VMD 


r 


BR th x 10 4 


BRe X p x 10 4 


77 — > 7T°77 (eV) 


1.24 x 1(T 3 


4.5 x 10~ 4 


0.26 


0.28 


2.1 


2.7 ±0.5 


r\' — > 7r°77(keV) 


7.7 x 1(T 5 


1.3 x 1(T 4 


1.29 


1.29 


65 


<8 (90% CL) 


r\' — > rjyy ( e V) 


1.4 x 1(T 2 


0.96 


48.8 


51.2 


2.6 





Table 1: Chiral-loop, LaM and VMD predictions for 77 — > n°yy, 77' — > n°yy and r\' — > T777. The total decay 
widths are calculated from the coherent sum of the LcrM and VMD contributions. The comparison between 
the predicted branching ratios and the present experimental values, if available, is also performed. 



In Ref. [11], the VMD prediction for 77 — > n°yy was calculated assuming equal p and ft) con- 
tributions and without including the decay widths in the propagators. In this case, these approxi- 
mations are valid since the phase space available prevents the vector mesons to resonate. However, 
for 77 ' — > 71° 77, the phase space allowed permits these vectors to be on-shell and the introduction 
of their decay widths is mandatory. For this reason, we include, for all the three cases, the decay 
widths in the vector meson propagators. 



4. LaM prediction 

An "all order" estimate of the scalar meson exchange effects to the processes under study can 
be achieved in the LaM where the complementarity between this model and ChPT can be used 
to include the scalar meson poles at the same time as keeping the correct low-energy behavior 
expected from chiral symmetry. This procedure was applied with success to the related V — > P°P Q y 
decays [17]. The t?o(980) enters into the calculation of r\ — > 7T°77and rj' — > n°yy, more intensively 



in the latter case on account of phase space, while the a(600) and /o(980) do the same in r\' — > 77 yy, 
although only the first contributes in a substantial way. Taking into account the scalar meson effects 
in an explicitly way does not provide a noticeable improvement with respect to the chiral-loop 
prediction, except for the case of 77' — > rjyy where the a contribution turns out to be considerable. 



However, the details of this calculation are involved and will be described elsewhere [18] 



5. Preliminary results 

The preliminary results of our analysis are shown in Table [j], where the predictions of chiral 
loops, the LaM, which replaces the former when scalar meson poles are incorporated, VMD, and 
the total decay width and branching ratio for the three processes are included. The comparison 
with the experimental results, if available, is also displayed. The total decay width is the result of 
adding the LaM and VMD contributions coherently. For the numerical results, we use f n = 92.2 
MeV, \g\ = 4.2 from the present value of p — > KK, and q> P = (40.4 ±0.6)° [Q for the 77-77' mixing 
angle. For 77 — > 7i°yy, our calculation agrees with the "all-order estimate" of Ref. [[ll]] and the 



more involved analysis of Refs. ]14j, |15J], thus giving support to our approach as a starting point 
for the determination of the other two processes. For 77' — > n () yy, the intermediate vector meson 
contributions dominate and the scalar meson effects are seen to be negligible. The ft) contribution 
prevails with a 80.2% of the total VMD signal, while the p contributes with a 4.6%. The predicted 
branching ratio appears to be one order of magnitude bigger than the old experimental upper bound. 
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Therefore, a new measurement would be welcome. Finally, for r\' — > f]yy, the VMD contribution 
also dominates but the scalar meson effects seem to be sizable, in particular those related with the 
a meson. The interference term is constructive. The p, ft) and <p contribute with a 59.9%, 15.8% 
and 1.6%, respectively, while the LaM calculation enhances by two orders of magnitude the chiral- 
loop prediction. Since G-parity does not apply to this case, the loop of charged kaons is suppressed 
and only the charged-pion loop plays a role. 

In summary, the calculated branching ratios for r\' — ^ n°yy and tj' — > r\ 77 give values which we 
consider are large enough to be measured in the near future by several experimental collaborations. 

Disclaimer 

For reasons of limited space, we do not include any of the two-photon invariant mass spectra 



which will be accessible in Ref. (18]. 
References 

[1] G. Amelino-Camelia, F. Archilli, D. Babusci, D. Badoni, G. Bencivenni, J. Bernabeu, 

R. A. Bertlmann and D. R. Boito et al, Eur. Phys. J. C 68, 619 (2010) [arXiv: 1003.3868 [hep-ex]]. 

[2] H. -B. Li, J. Phys. G G 36, 085009 (2009) [arXiv:0902.3032 [hep-ex]]. 

[3] S. Prakhov, B. M. K. Nefkens, C. E. Allgower, V. Bekrenev, W. J. Briscoe, J. R. Comfort, K. Craig 
and D. Grosnick et al, Phys. Rev. C 78, 015206 (2008). 

[4] D. Aide et al. [Serpukhov-Brussels-Annecy(LAPP) and Soviet-CERN Collaborations], Z. Phys. C 25, 
225 (1984) [Yad. Fiz. 40, 1447 (1984)]. 

[5] S. Prakhov et al, Phys. Rev. C 72, 025201 (2005). 

[6] K. Nakamura et al [Particle Data Group Collaboration], J. Phys. G G 37, 075021 (2010). 

[7] M. Unverzagt [Crystal Ball at MAMI Collaboration], Nucl. Phys. Proc. Suppl. 198, 174 (2010) 
[arXiv:0910.1331 [hep-ex]]. 

[8] P. Gauzzi [KLOE-2 Collaboration], J. Phys. Conf. Ser. 349, 012002 (2012). 

[9] B. Di Micco et al. [KLOE Collaboration], Acta Phys. Slov. 56, 403 (2006). 

[10] D. Aide et al. [Serpukhov-Brussels-Los Alamos-Annecy(LAPP) Collaboration], Z. Phys. C 36, 603 
(1987). 

[11] L. Ametller, J. Bijnens, A. Bramon and F. Cornet, Phys. Lett. B 276, 185 (1992). 
[12] J. N. Ng and D. J. Peters, Phys. Rev. D 47, 4939 (1993). 

[13] Y. Nemoto, M. Oka and M. Takizawa, Phys. Rev. D 54, 6777 (1996) [hep-ph/9602253]. 

[14] E. Oset, J. R. Pelaez and L. Roca, Phys. Rev. D 67, 073013 (2003) [hep-ph/02 10282]. 

[15] E. Oset, J. R. Pelaez and L. Roca, Phys. Rev. D 77, 073001 (2008) [arXiv:0801.2633 [hep-ph]]. 

[16] R. Escribano, P. Masjuan and J. J. Sanz-Cillero, JHEP 1105, 094 (2011) [arXiv: 101 1.5884 [hep-ph]]. 

[17] R. Escribano, Phys. Rev. D 74, 114020 (2006) [hep-ph/0606314]. 

[18] R. Escribano and R. Jora, work in preparation. 

[19] F. Ambrosino, A. Antonelli, M. Antonelli, F. Archilli, P. Beltrame, G. Bencivenni, S. Bertolucci and 
C. Bini et al, JHEP 0907, 105 (2009) [arXiv:0906.3819 [hep-ph]]. 



6 



